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Dynamical systems driven by a general Levy stable noise are considered. The inertia is included 
and the noise, represented by a generalised Ornstein-Uhlenbeck process, has a finite relaxation time. 
A general linear problem (the additive noise) is solved: the resulting distribution converges with time 
to the distribution for the white-noise, massless case. Moreover, a multiplicative noise is discussed. 
It can make the distribution steeper and the variance, which is finite, depends sublinearly on time 
(subdiffusion). For a small mass, a white- noise limit corresponds to the Stratonovich interpretation. 
| On the other hand, the distribution tails agree with the Ito interpretation if the inertia is very large. 

■ An escape time from the potential well is calculated. 
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I. INTRODUCTION 

A prominent feature of the stable Levy processes is the existence of algebraic, long tails of the probability distri- 
butions of the form |x| -a_1 , where < a < 2 is a stability index. As a consequence, the moments, in particular the 
variance, are divergent. The diffusion process in such systems is called an "accelerated diffusion" and the relative 
transport rate may be described by a time-dependence of the fractional moments instead of the variance. However, 
, 1 * . there are indications that in some physical problems the distribution tails fall faster than for the pure Levy flight. 
In the field of the economic research, the indexes 2.5 - 4 were observed in financial data [l|; it has been suggested 
that such values of the index arise when the trading behaviour is performed in an optimal way 0|. The probability 
-4_J ■ distributions of the hydraulic conductivity in the porous media seem to obey the power law with the index 3.5, while 
the atmospheric turbulence studies yield even larger index for the wind field Q. Such slowly falling algebraic tails 
are predicted by the Langevin equation with the Levy stable noise - in a sense of the stationary solution - when 
one introduces an appropriate deterministic potential. It has been demonstrated by Chechkin et al. [3] that the 
stationary distribution tails of the form j;-"- 2 ™- 1 result from the potential ~ x 2m+2 . Also temporal characteristics 
of the system may influence the asymptotic shape of the distribution. This happens if, for a jumping process, long 
jumps are penalised by a short waiting time. The finite variance is observed for such jumping processes as the Levy 
1— ~ 1 walk 5] and the kangaroo process with a Levy distributed jumping size 6]. 

, The Levy stable processes are often connected with complex phenomena for which the power-law shape of the 
distribution tails [7| is typical, as well as a complicated structure of the medium. It is the case for the porous media, 
V/~j | plasmas and fractal (multifractal) structures [3j, |8[ . Therefore, a nonhomogeneity must often be taken into account in a 
I/"") ■ dynamical description, both as a deterministic potential and as a multiplicative noise. Descriptions of the diffusion on 
C^ - 1 fractals involve the variable, power-law diffusion coefficient [1, |T(| ■ Also the other topologically complicated systems 
with long jumps, the folded polymers, require a variable diffusion coefficient to describe the transport fllj . Moreover, 
formalisms with the multi plic ative Levy noise can describe the second order phase transitions [12j and the dynamics 
of two competing species (13(. A nonlinearity of the Langevin equation makes the stochastic process different from 
the pure Levy motion. In particular, variance may be finite for a system driven by the multiplicative Levy noise [l4j . 
Generally, the variance rises not only linearly with time but also faster or slower than that, i.e. the diffusion may be 
anomalous. In the case of Ref.[14|, motion is subdiffusive. The above approach includes the white noise. However, 
a Markovian description of a realistic system is an idealisation, valid only if the time scale of fast variables is short, 
compared to the time scale of the process variable. A procedure of the fast variables elimination produces correlations: 
they are present even if the original system is Markovian 15]. It has been demonstrated for the Gaussian noise that 
characteristic time scales of the fast variables are important even if the variables themselves are eliminated [l6j : this 
finding suggests using a coloured noise in a stochastic description rather than the white noise. Effects related to the 
correlations are important, for example, for such problems as fluctuations of a dye laser light [l7j and a narrowing 
of the magnetic resonance lines (l8| . Importance of the finite correlation time for noise- induced phase transitions 
was emphasised in Ref. 19]; an increase of that time favours disorder and prevents the formation of an ordered state. 
Introducing the white noise as a limit of the finite correlation time means that the stochastic integral should be 
interpreted in a Stratonovich sense (20j . On the other hand, effects of the finite inertia should be taken into account. 
If the relaxation time associated with the inertia is large compared to the correlation time, the Ito interpretation 
comes into play [2lJ. That effect of the inertia, opposite to the correlations, was demonstrated in Ref. (22| : it modifies 
the front propagation by suppressing the external multiplicative, white noise influence on the velocity of fronts. 
The Ito-Stratonovich dilemma becomes especially interesting for a < 2 since then - when we consider the white 
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Levy noise and neglect the inertia - the very existence of the variance depends on the particular interpretation of the 
stochastic integral. This problem is important for the diffusion since the infinite variance, which means the infinite 
propagation speed, is unphysical in most cases. How do the finite noise relaxation time and the inertia modify slope of 
the distribution? We address this question in the present paper and discuss consequences for the diffusion. In Sec. II a 
linear problem involving the additive noise is considered. Sec. Ill is devoted to the multiplicative noise; both limiting, 
analytically solvable cases and numerical solutions are discussed. Moreover, the escape time from a potential well is 
calculated. Results are summarised in Sec. IV. 



II. ADDITIVE NOISE 



We consider a linear problem which is defined by the following system of the Langevin equations for x, v and £: 



mv(t) = f Q - 0v(t) - Xx + 7 f(t) 
x(t) = v(t) 

d£(t) = -i£(t)dt + dL(t), 



(1) 



where is a damping coefficient. The stochastic force, dL(i), is the symmetric and stable Levy process, characterised 
by the stability index a G (0,2]. Special cases of the system ([1]) were considered by several authors. The velocity 
distribution for the white noise without a potential was obtained in Ref.(23|, the linear force case was discussed in 
Ref.[24| and the white noise case with the inertia for a = 1 in Ref.[25[. Moreover, the asymmetric Levy distribution 
was introduced in Ref.[26]. On the other hand, stochastic collision models may lead to the Levy statistics. The 
Fokkcr-Planck equation with the additive noise predicts, in the limit of small mass, an equilibrium in the form of the 
Levy distribution [13] . The case a = 2 corresponds to the normal distribution. Then the third equation ([T]) describes 
the standard Ornstein-Uhlenbeck process with the covariance 
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therefore 7 determines a correlation time, 1/7. A generalisation of the Ornstein-Uhlenbeck process for a < 2 implies 
an infinite covariance for any time. However, the parameter 7 can still estimate the noise relaxation time. One can 
modify the covariance definition [2^, to get a convergent quantity which behaves with time similar to Eq.©. 
On the other hand, the covariance becomes finite when one introduces a truncation of the Levy distribution [30j. 
Properties of such a dynamical system are similar to the system without the truncation for an arbitrarily large time 
[3ll | and the parameter 7 measures the correlation time. Values of the process dL(t) are given by the characteristic 
function p(k) — exp(— K a \k\ a ) (K > 0). The fractional Fokker-Planck equation 
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determines the probability density distribution p(x, v, £; t\xo, vq, £0; 0) and the fractional Weyl derivative is defined 
by its Fourier transform, J~[ a ^ x \ a f{x)} — —\k\ a f{k). We will evaluate the density of x, p(x,t), directly from the 
stochastic equation (JTJ. We restrict our analysis to the case of a relatively weak potential; more precisely: let 
/3 2 /m 2 — 4A/m = A 2 > 0. First, we need to evaluate the stochastic trajectory x(t). The solution of Eq.([T]) produces 
the result 



:(t) = A+ f ' h(f)L(t')df, 
Jq 



(4) 



where 



^e-^ 2m sinh-t+^ 
A 2 mA 



2A 6XP ^ 2m 1 /3/m + A (3/m - A 



(5) 



and 
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m 7 2 — 10 + X 



(6) 
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Two simple special cases are distinguished. In the absence of inertia (m = 0), we have an adiabatic problem of the 
particle subjected to the linear force and the coloured Levy noise. Then Eq. (JXJ) yields 
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Secondly, for the case of a free-particle (A = 0) we obtain 



j-P/m (3* C 



7 (t-t') 



1 



(t-t')/m 



(J) 



(8) 



(9) 



Eq. (|9]) implies the Levy distribution with the same stability index as the driving noise L and a translation parameter 
which coincides with fx and for large time equals either /o/A (A ^ 0) or fo(t — m/ (3) (A = 0). Since the dependence 
of the density distribution on /o is trivial, we assume in the following fo — and vq = 0. Then the distribution is 
symmetric for any time. The inverse Fourier transform can be conveniently expressed in a form of the Fox function 
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The characteristic function of p(x, t) directly follows from Eq.© [H, [j^]: 

p(k,t) = (c lkx ^) = c lkfl (exp(ifc / f 2 {t')dt')) = c lkfl exp ( -K a \k\ a I \f 2 (t')\ a dt' 



p(x,t)=Nf(t)H. 



2,2 



f(t)\x\ 



(l-l/a,l/a),(l/2,l/2) 
(0,1), (1/2, 1/2) 



(10) 



where f(t) — K a f Q \ f2(t')\ a dt' . By introducing a new variable r = t — t' we have f(t) — K a (f + f^) and the 
exponentials in the second integral can be dropped for any T>1. Therefore f(t) = t//3+const for large times. Since 
p{x,t) ~ |a;| -1 ~ a for \x\ — > oo, variance and all higher moments of the distribution (I10[) are divergent, as well as the 
average if a < 1. A relative expansion rate can be quantified by the fractional moments of the order S < a, (|a;|*)(t); 
they are given by the Mellin transform from the Fox function [5( . The final expression reads 



(\x\*\- If (*)*'* r M/") r ( 1 + ^ 
{l 1 ^ a J{> T(-S/2)T(l + 8/2)- 



(11) 



Consequently, in the limit of large t the fractional moments decrease with the damping coefficient /3 and the expansion 
rate is large for small a. 

As an example, let us consider the case a = 1 for which results take the transparent form. This particular value of 
the stability parameter corresponds to the well-known Cauchy distribution; it was considered in Ref. [25| for the white 
noise. For 7 > /3/m, a straightforward calculation yields the expression for the apparent width of the distribution 
p(x,t): 
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In the limit of the large time, inertia and noise relaxation time are responsible for a time-shift which is negative and 
rises with m and I/7. 



III. MULTIPLICATIVE NOISE 



By introducing a multiplicative noise we take into account that the influence of the random component of the 
dynamics depends on the dynamics itself. The one-dimensional case is given by the Langevin equation 

mv(t) = -dV(x)/dx - Pv(t) + 7 G(x)£(f) 
x{t) = v(t) 

d£,{t) = -rt(t)dt + dL(t). (13) 
In the following we assume the noise intensity in the algebraic form, G(x) = \x\~ e / a (a + 9 > 0). 
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The simplest problem involves the white noise (7 — > 00) and neglects the inertia. The latter condition means that 
mass is small compared to the damping parameter /? and the strength of the noise. The case of the normal distribution, 
a = 2, is well-known (35l - [37j . A stochastic integral in the Langevin equation is not completely determined for the 
uncorrelated noise, since it is not clear whether the dynamical variable in the function G(x) should be evaluated at 
the time before the noise acts, after that, or somewhere in between. Two interpretations are of particular importance. 
In the Ito interpretation, the i's component of the discretized stochastic integral is G[x(ti-i)][£(ti) — £(t*-i)], whereas 
the Stratonovich interpretation includes both the beginning and the end of each interval: G{[x(ti-i)+x(U)]/2}[£(ti) — 
£{U-i)]- Both assumptions result in a different Fokker-Planck equation but a difference resolves itself solely to a drift 
term (the spurious drift) which can be eliminated by an appropriate modification of the deterministic potential. For 
this reason, a physical relevance of the Ito- Stratonovich dilemma is disputed [36]. Nevertheless, physical implications 
of both interpretations may be different. For example, phase transitions due to instabilities of the disordered phase in 
the framework of the Ginzburg-Landau model take place only for the Stratonovich interpretation 38]. Rules of the 
ordinary calculus are valid in the Stratonovich formalism, in contrast to the Ito interpretation. 

The meaning of both interpretations becomes more transparent when we take into account the finite correlations 
and inertia. First of all, the memory effects favour the Stratonovich interpretation since it constitutes the white- 
noise limit of the correlated processes [2(| ■ Inertia acts in the opposite direction. It was demonstrated [2l[ - by the 
estimation of the velocity moments and using the Ito formula - that if the inertia relaxation time goes to zero faster 
than the noise correlation time, the Stratonovich interpretation is valid. The opposite limit produces the Ito result. 
If both time scales are comparable, neither of the above interpretations is valid. We will demonstrate that similar 
conclusions can be drawn for the general Levy stable processes. However, for a < 2 methods of Ref. 21] cannot be 
applied since the moments are divergent and the Ito formula is unknown. 

We take into account only Ito and Stratonovich interpretations of the stochastic integral. The Ito interpretation 
applies, beside the systems with large mass, to discrete problems; it is commonly used in the perturbation theory 
[391 ]. However, there are indications that other interpretations are also important. For example, it was recently 
experimentally demonstrated that description of the Brownian motion in the presence of gravitational and electrostatic 
forces requires a backward integral (anti-Ito interpretation) [40| . 

A. White-noise case without the inertia 

To study a diffusion process, we consider a free particle, V(x) = 0. In the limit m — >• and 7 — > 00, Eq. (TT3^) 
becomes a single Langevin equation of the first order: 

dx{t) = \x{t)\- e/a dL(t), (14) 

where, for simplicity, we assumed = 1. In the Ito interpretation it corresponds to the Fokker-Planck equation 0] 

^ Pl ( X) t) = K a -^[\x\- e Pl (x,t)}, (15) 

which differs from the equation for the additive noise by an algebraic term under the fractional derivative. This 
particular form of the multiplicative factor suggests simple scaling properties and a possible similarity of the solution 
to Eq. fTU)) . Indeed, an asymptotic solution of Eq. fPo)) can be found in the form a(t)H 2 , 2 [a(t)x]. The procedure is 
the following (4l| . First, we insert the above expression to the Langevin equation and take the Fourier transform, 
which also has a form of the Fox function but of a higher order. Expansion of the Fox functions in powers of k and 
neglecting the terms of the order |fc| 2Q+ and higher yields a simple differential equation for the function a(t) and 
allows us to determine some Fox function coefficients. Finally, we obtain the solution 



pi(x,t) =Na(t)H 2 1 < 2 1 



a(t)\x\ 



i 1 ~ STe)' ("2,^2) 
(h,Bx), (1 - 5^3, 2^5) 



(16) 



where a(t) ~ t 1 /( a + e ) and the coefficients (02, A2) and (61, B\) are arbitrary. The asymptotic form of the solution is 
the same as for the driving noise, 

Pl {x,t)r.\ X \- 1 - a . (17) 

Therefore, for the Ito interpretation the variance is always divergent which implies accelerated diffusion. Fractional 
moments can be evaluated similarly to the additive noise case; a straightforward calculation yields 

{\x\ s } ~t s ^ a+e \ (18) 
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FIG. 1: (Colour online) Probability density distributions at t = 1 for a — 1.5, 8 = 2, 7 = 100 and /3 = 1. The curves correspond 
to the following values of m (from left to right): 5 x 10 5 , 5 x 10 4 , 10 4 , 10 3 , 0.01, 0.1, 10, 1. The red dashed line marks the 
depedence x~ 4 ' 5 and the green dashed line (at the left side) the dependence a; -2 ' 5 . Each curve was obtained by averaging over 
10 8 trajectories. 



The multiplicative noise parameter 9 can both strengthen [9 < 0) and weaken (9 > 0) the time-dependence of (\x\ s ), 
compared to the case of the additive noise. The undetermined coefficients do not influence the functional dependences. 
The same method of solution can be applied in the presence of the linear deterministic force F(x) = —Xx 

Results for the Stratonovich interpretation are qualitatively different since the decline of the noise intensity with 
x may compensate the effect of the long jumps. The technical advantage of this interpretation, for one-dimensional 
systems, consists in a possibility of applying rules of the ordinary calculus. This property is strict for a = 2 |39| . 
In the general case, the noise distribution must be truncated, a requirement that is obvious for the linear systems 
[42j . However, it was numerically demonstrated that in practice cases with the distribution without any cut-off also 
comply with rules of the ordinary calculus if the system is nonlinear [bl . |42| . Then we may define a new variable, 

y(x)= K{ " +e) \x\ 1+(> /*sgn(x), (19) 

which transforms Eq. (ll4[) to the equation with the additive noise. The asymptotic form of the solution [lij 

p s (x,t) ^ t a / {a+e) \x\- 1 - a - 9 (M^oo) (20) 

implies that variance may be convergent. It takes the form (a; 2 ) ~ i 2 /( Q + ) on the condition a + 9 > 2. Therefore, 
diffusion is either anomalously weak - if the above condition is satisfied - or accelerated [43j . 



B. General case 



First let us consider the overdamped limit (the adiabatic approximation) by putting m = in Eq.(|13j). Equations 
take the form 

m - j\x\- e/a m 

d£(t) = -j£(t)dt + dL(t). (21) 
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FIG. 2: (Colour online) Same as Fig.l but for 7 = 1. The following values of m are presented (from left to right): 5 x 10 , 
5 x 10 4 , 10 4 , 10 3 , 10 2 , 0.01, 0.1, 1. 



After transformation of the process variable according to Eq. (fT9|) , we obtain a linear equation with the additive noise. 
Its solution reads y(t) — J Q f2(t')L(t')dt' , where fait') = (1 — e~ 7 ^ t_t ')//?, and the density distribution of y has the 
Levy form, Eq. fTU)) . Transformation to the variable x yields the final result: 



p(x,t) 



a+e/a 



K{l + d/a)f{t) 1 / a 



(1,1/a), (1,1/2) 
(1,1), (1,1/2) 



(22) 



where f{t) — K a J"' \f2(t')\ a dt' . The expansion of Eq. ((22|) in the fractional powers of \x\ 1 yields an approximation 
of the solution for large \x\ and the first term is of the form 

p(x,t) ~ f{t) a/{a+9) \x\~ 1 ~ a ~ e . (23) 

If a + 8 > 2, the variance is convergent and it can be exactly evaluated by using properties of the Fox functions, in 
particular, an expression for the Mellin transform. A straightforward calculation yields 
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fit) 



2 /(a+9) 



(24) 



Eq. (|22[) converges to the white- noise solution in the Stratonovich interpretation for any noise relaxation parameter 7 
if time is large or, for any time, if 7 -+ 00. 



Distributions for arbitrary 7 and m have been obtained by a numerical integration of the stochastic equations, 
Eq.(fT3"|). For that purpose, a second order difference approximation, called a Stormer method [44j, was applied to 
the first two equations. Since the resulting difference equations are implicit, the parabolic interpolation scheme was 
applied at each step (45j. The third equation was integrated by an Euler method and the noise term in the i-step 
was represented by r 1 / Q Li, where r was a time step 46]. We will demonstrate how the asymptotic shape of the 
distribution, for a given time, depends on m and 7. Another quantity of interest is a time dependence of the variance. 
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FIG. 3: (Colour online) Slopes of the tails, |a;| _M , of p(x,t) at t = 1 for a = 1.5, = 2 and (3 = 1, as a function of m. Two 
cases are presented: 7 = 100 (points) and 7 = 1 (squares). The limiting values, [is and fii, are marked by the horizontal lines. 
Inset: \i as a function of 7 for m = 5x 10 . 



Fig.l presents the probability density distributions as a function of the particle mass in the limit of the white noise 
at t = 1. The distributions widen with m up to m = 1 but then the trend goes into reverse. For the large mass the 
distributions have a form of the delta function accompanied by a little tail. The tails are algebraic, ~ and the 

slope fi diminishes with the mass. Two limiting values, [is = a + 9 + 1 and [ii — a + l, correspond to the Stratonovich, 
Eq. ([20| . and Ito, Eq. (flTj) . interpretations, respectively. The distributions for the case of the finite noise relaxation 
time (7 = 1) are presented in Fig. 2. They are similar to those for the white noise but the limiting slopes are not yet 
reached &tt = l. Slopes for all cases are put together in Fig. 3. The dependence n(m) is flat for 7 = 1 whereas for the 
white- noise case large values of fx dominate and there is a rapid transition to fii. Variance is finite (/1 > 3) except for 
the very large m. This case is separately presented in Fig. 3: [i rises with the noise relaxation time from the Ito value 
for the white noise, [i = 2.5, to /i = 3, where it saturates. 



Diffusion properties of the system are determined by a long-time behaviour of the variance. The case corresponding 
to the short noise relaxation time is presented in Fig. 4. If m is very small, the variance assumes the form (a; 2 ) ~ t 2 l 
for the large time. The slope becomes slightly larger if m is not infinitesimal; it equals 0.68 for all m > 0.1. 
Therefore, all cases indicate a sub-linear time dependence for the large time: the diffusion process is anomalously 
weak (subdiffusion) . On the other hand, if time is not very large, (x 2 )(t) exhibits a plateau which widens with m. 
Moreover, the curves reveal a stepwise pattern which can be attributed to a competition between the expansion and 
the attraction to the origin. Such a behaviour of the curves in Fig. 4 is a clear consequence of the lack of memory. 
For 7=1 the dependence (x 2 )(t) - presented in Fig. 5 - is smooth; it assumes the asymptotic shape i 068 for large m, 
similar to the previous case. If m is close to zero, variance is given by Eq. (|2"4"]) . 

Also properties of more complicated systems are modified when we take into account the finite relaxation time and 
inertia. Let us consider the following potential 



V(x) = ^ - fx 2 



(25) 
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FIG. 5: (Colour online) The same as Fig. 4 but for 7=1. The curves correspond to the following values of m (from bottom to 
top on the right side): 0.01, 0.1, 1, 2, 5, 10, 15, 20. 
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FIG. 6: (Colour online) MFPT as a function of 9 for finite m and 7. The black solid lines correspond to the following values 
of m (from bottom to top): 0.001, 0.01, 0.1, 1, 10 and 7 = 100. Dependence on 7 for m = is marked by green dotted lines: 
7 = 1000, 10, 1, 0.5, 0.2 (from bottom to top). Other parameters: a — 1.5, ft = 1, A = 1 and B = 0.1. Result for the 
Stratonovich interpretation is marked by the red solid line without symbols. 

which has the double well shape. The mean first passage time (MFPT) is a quantity of particular importance [47j : 
it was studied in the context of the Levy stable processes in Ref. (48l . |49| . The case of the multiplicative noise was 
discussed in Ref.[i3|; it was demonstrated that MFPT decreases with 9 but the rate depends on the particular 
interpretation of the stochastic integral. In this paper we calculate MFPT for finite m and 7 by integration of Eq. (fT5)) 
with the absorbing barriers at x(0) = —^JB/A and x = 0. The latter boundary condition is nonlocal due to the 
jumps (49[. The resulting MFPT, as a function of 9, is presented in Fig. 6. All the curves fall since the effective depth 
of the potential decreases with 9. MFPT rises with m, because of the increasing attraction to the origin, and becomes 
flat. A similar effect is observed for the decreasing 7 (stronger memory) since then the intensity of the driving noise 
is smaller. In the white noise limit, 7 — > 00, the Stratonovich result is recovered. 



IV. SUMMARY AND CONCLUSIONS 

Wc have studied a one-dimensional dynamics of a massive particle subjected to the general Levy stable noise, both 
additive and multiplicative. The driving noise has been represented by the generalised Ornstein-Uhlenbeck process 
and then the finite noise relaxation time has been taken into account. In the linear case, the dynamical variable x is 
governed by the Levy distribution and the parameter a is the same as for the driving noise. Therefore, diffusion is 
always accelerated. Distribution converges with time to the white-noise and massless case. Fractional moments rise 
with time; the rate decreases with the stability index a and the damping coefficient j3. Inertia and noise relaxation 
time influence the rate of convergence to the asymptotic distribution. 

The ^-dependence of the multiplicative noise modifies the distribution. Slopes of the tail depend on the multi- 
plicative factor G(x), which was assumed in the algebraic form, and variance is finite if G(x) falls sufficiently fast. 
Variance rises sub-linearly with time for t ^S> 1 which indicates the subdiffusion. Those conclusions are valid for any 
noise relaxation parameter 7. The limit 7 — > 00 is of particular importance; the distribution in this limit coincides 
with Eq. (|20[) . Therefore, the limit of the Langevin equation driven by the generalised Ornstein-Uhlenbeck process 
produces the same result as the formal variable change in the Langevin equation for the white-noise case. The influ- 
ence of inertia is more subtle. It favours an expansion of the distribution if m is small but for large m distribution 
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shrinks to the delta function. However, even in the limit m — > oo a little tail remains and it makes the variance 
divergent. In the white-noise limit, that tail agrees with the distribution in the Ito interpretation. On the other hand, 
the Stratonovich interpretation is valid for the small mass. Those conclusions are similar to the case of the normal 



distribution [21|. Since slowly falling tails have been encountered only for the extremely large masses, convergent 
variance is by no means exceptional for the Levy stable processes: it emerges if intensity of the multiplicative noise 
diminishes sufficiently fast. The finite noise relaxation time and inertia affect the barrier penetration: the calculated 
MFPT rises with both the memory parameter I/7 and the particle mass. In the white-noise limit MFPT converges 
to the Stratonovich result. 

The above analysis demonstrates that the Langevin formalism with the multiplicative Levy noise predicts heavy, 
algebraic tails of the probability density distribution and the index /x can assume arbitrarily large values. As a 
consequence, moments of an arbitrarily high order may be convergent, fi depends not only on a and 9, as it is the 
case for the massless particle, but also on the inertia. Those conclusions suggest that the presented formalism may 
be well suited to describe processes characterised by a variety of the algebraic slopes of the distribution [l|-[3| . In the 
field of finance, a traditional Black-Scholes model of option pricing, which includes the additive Gaussian noise, can 
be generalised by introducing the Levy flights. Need of such a generalisation is obvious [5(3] but, since variance of the 
additive Levy process is infinite, a truncation of the distribution becomes necessary. On the other hand, the first-order 
equation, like the Black-Scholes equation, with the multiplicative noise predicts sufficiently steep distribution slopes 
to ensure the finite variance also without any truncation if the stochastic integral is understood in the Stratonovich 
sense. The present paper justifies this interpretation for the first-order stochastic equations: it demonstrates that 
distribution slopes are robust in respect to the noise relaxation time - which is always finite for realistic problems - 
and the white-noise limit exists. 
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